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The general one-dimensional "log-sine" gas is defined by restricting the positive and 
negative charges of a two-dimensional Coulomb gas to live on a circle. Depending on charge 
constraints, this problem is equivalent to different boundary field theories. 

We study the electrically neutral case, which is equivalent to a two-dimensional free 
boson with an impurity cosine potential. We use two different methods: a perturbative 
one based on Jack symmetric functions, and a non-perturbative one based on the thermo- 
dynamic Bethe ansatz and functional relations. The first method allows us to compute 
explicitly all coefficients in the virial expansion of the free energy and the experimentally- 
measurable conductance. Some results for correlation functions are also presented. The 
second method provides in particular a surprising fiuctuation-dissipation relation between 
the conductance and the free energy. 
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1. Introduction 

The general 2D classical Coulomb gas with charges restricted to live on a circle is a 
recurrent problem in several areas of theoretical physics. These include random matrix 
theory [|l|, impurity problems (like the Kondo effect 0, resonant tunneling in quantum 
wires [|] and between quantum Hall edge states and dissipative quantum mechanics 

H,[7|,|8|j9|,|T0| ■ In this Id "log-sine" gas, the charges interact with a long-range interaction 
proportional to the log of the sine of the separation. 

Two particular cases of this Coulomb gas have been solved analytically. The gas 
with only one type of charge (the Dyson gas) is related to eigenvalue statistics for circular 
ensembles [Q , and can be addressed by elementary methods. When there are two types of 
charges required to alternate in space [H , the gas is related to the Kondo problem , and 



therefore can be addressed indirectly by the Bethe ansatz solution of the latter []Tl |. 

We present in this paper two methods to address the general model with two types 
of charges. The first method is direct, and uses the recently-developed technical tool of 



Jack symmetric functions |[T^,|T3|,|IJ|T^] . These have been used extensively in recent works 
on Green functions for the Calogero- Sutherland model [p!6| , p!7|JT^ , |T9| , |20|j2l]] . The second 
method is indirect and uses the solution of the boundary sine-Gordon theory via exact 
S matrices and the thermodynamic Bethe ansatz (TBA) [p3 ,24|. 

Technically, the combination of the two methods allows us to compute a rather large 
number of quantities, including dynamical properties. When they overlap, they can be 
compared, leading to interesting relations between two different active areas of mathe- 
matical physics: models and the TBA. Physically, our solution can be applied to a 
range of interesting problems, including the experimentally-measurable resonant tunneling 
between quantum Hall edge states (which was derived using the TBA in 0), and the case 
of dissipative quantum mechanics which will be discussed elsewhere. 

We start by describing this model as a l-|-l-dimensional field theory with an impurity. 
Consider a Gaussian model on an infinite cylinder with action 

-| POO pR pR 

A= — / / da dr [{dr(l)f + {d^cpf] + 2g dr cos [(3 (p{a = 0,t)] . (1.1) 
2 J-oo Jo Jo 

With periodic boundary conditions in the r-direction, this is equivalent to a one- 
dimensional quantum problem at non-zero temperature 1/R, with an impurity at a=0. 
The impurity coupling g has a dimension, so this problem is not conformally invariant and 
the interaction induces a flow from the free boson on an infinite cylinder in the UV {g=0) 
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to two decoupled free bosons on two half-cylinders with Dirichlet boundary conditions 
at their boundary in the IR {g large). A convenient quantity describing this flow is the 



"g- factor" discussed in |j2^ , whose logarithm is equal to the impurity entropy (the contri- 
bution to the entropy which is independent of the length of the cylinder). This "g- factor" 
(which we prefer to denote u here) is w = 1 in the UV and uj = t~^l'^ in the IR, where we 
deflne 

t 



47r 



We can study (|1 . 1|) by naive perturbation theory, which exhibits the relation with a 
Coulomb gas. Deflning as usual the partition function as Z = /[d^] we introduce 

^_ Z{g) 
- Z{g = 0)- 

Using the free fleld propagator on an inflnite cylinder 
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(where k is a renormalization constant) we obtain by the standard perturbation expansion 
in powers of g 
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The sum (|1.2| ) is the grand canonical partition function for a classical two-dimensional 
Coulomb gas with two species of particles (with positive and negative charges) that lie on 
a circle of radius R, and which is electrically neutral. The parameter g is the fugacity of 
charges, while (3 controls the combination charge ^ x inverse temperature. In sect. 2 we 
will show how these integrals can be explicitly evaluated by expanding the integrands in 
terms of Jack polynomials. 

We can also reformulate this as a boundary problem (i.e. on the half- cylinder) following 
2^ . We introduce two new flelds 



a, r) = — = [0((j, r) + ^(-a, r)] 
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so the action reads now 



-I poo pn 

A=-- / dadr [{dM'' + {dM'' + (5r<^o)' + {d^''] 
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(1-3) 



+2g / dr cos 







One can of course obtain the perturbative expansion ( p. .21) from ( p. .31) by using the propa- 
gator on the half-cyhnder. In the UV both fields have Neumann boundary conditions and 

1/2 

uj = (I) while in the IR the odd field has still Neumann boundary conditions while 
the even has Dirichlet boundary conditions and = Notice that in the reformulation 
the absolute values of to have changed; this is presumably due to Jacobian terms in the 
functional integral when redefining the new fields. However, the ratio ^^^^ = — = 

° o ' UJIR UN 

remains a constant. 

Using this reformulation, the free energy JF = —Tin 2 can be obtained non- 



perturbatively using the thermodynamic Bethe ansatz [p3| , p^ . The corresponding analysis 
appears in ||2^ for t integer. We will discuss this in sect. 4, where we will also derive a set 
of functional equations which Z satisfies. The results of TBA depend on a dimensionless 
variable T/Tb, where 

f-i 

g = k'Tj^' , 

and is an unknown renormalization constant. The non-perturbative free energy also 
contains a non-analytic term T independent of T, which the detailed analysis gives as 

y= (14) 

2cos(7r/2(t- 1)) ^ ' 

There is also a shift because we have T{g = 0) = 0, whereas the TBA is defined so that 
^TBA{g oo) — > 0. Thus 

Ttba = -T\nVt + y. (1.5) 

This allows us to obtain the exact behavior of the free energy at large g (large Tg ) , because 
the the power series must be precisely balanced by JF, and thus 
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Z exp ( T—]^/ TTT ) (1-6) 



V2cos(7r/2(t-l)) 

in this limit. We will see in sect. 2 that this behavior also follows from the expansion ( |1.2| ). 
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This two-dimensional local field theory problem can also be reformulated as a one- 
dimensional non-local field theory on a circle of circumference R by integrating out the 
"bulk" degrees of freedom. Let us consider for instance the second (boundary) point of 
view and forget about the odd field which totally decouples. After integration the even 
action reads 



j^bdry 



TT 



„ , / ardT 7T + 2g dr cos 

'^R^ Jo Jo [sin f(r-r')] 
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We therefore have a one- dimensional model with a sine- Gordon- type interaction, where 
the Gaussian part has a l/r^ interaction. Related models have been considered in P^,ES 



2. Exact solution of the Coulomb gas problem using Jack polynomials 

In this section we derive all the coefficients of the perturbative expansion of the parti- 



tion function. The result (|1.2|) can be easily recast after a change of integration variables 
into 

oo 



(2.1) 
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where we have set 



and 
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dz, dz\ \ (A(2)A(2))V*(A(2')A(^'))iA 



2mz, 2mz\) [1, - - z'^z,)^!^ 



where A(2;) is the n-variable Vandermonde determinant A(2) = ni<j(-2j ^ -^i)? z^ 



To evaluate this integral, we expand the integrand in terms of Jack polynomials 
n 7137^ = E^^(«) Pxir,a)P,is,a). 
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2.3) 



The function P\{r,a) is a symmetric polynomial in the set of variables (ri, r2, . . . , r„) 
which depends on a rational number a. The subscript A is a partition of an integer; this is 
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conveniently labeled by a Young tableau; e.g. the partition 5 = 2 + 2 + 1 is labeled by the 
tableau with two boxes in the first row, two boxes in the second, and one in the third. The 
polynomials P\{x, a) vanish if the number of parts /(A) of the partition A (i.e. the number 
of rows of the tableau) is greater than the number n of variables, so the sum runs over all 
partitions of all integers with /(A) < n. The Jack polynomials have the useful property 
that their orthogonality relation involves the Vandermonde determinant: 



n 



2i7{-z,, 



iA{z)A{z)r Px{z, a)P^{z, a) = 6x,^Nx{a). 



(2.4) 



Hence the value of the integral follows (here a 

1 



1/t) : 



^2n 
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The numerical coefficients in the foregoing expression are expressed as a product over the 
boxes of the Young tableau associated with the partition A [|4|,|12|,|T^ . We have 

'j-l + i(n-z + l 



n 



j + jin-i) 



(2.6) 



where 



r(in+l) 



and s = is the box of the tableau at the intersection of the j*^ column and i^^ line. 

One can write this product compactly using gamma functions. Two convenient ex- 
pressions of follow, depending on whether one uses partitions or their conjugates (given 
by interchanging the rows and columns of the Young tableau). One obtains 



^2n 



ni 
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and alternatively, using the conjugates. 
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where Ai+i < A^. 
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Consider for instance the case n = 1 : (|2.7| ) reads then 

2 



The sum converges only for t > 2, so that the UV dimension of the perturbing operator 



X 



^ < \ BlHl- This coincides with the domain where the integrals in ( |1 . 2| ) are UV 
convergent. (We have no problem with IR divergences because we are on a circle.) The 
sum in ( |2.9| ) can be done explicitly; one can also obtain its value by a direct treatment of 
the integral avoiding Jack functions. The result is 

Z,(t) = i^i^. (2.10) 

r2(i-i/t) ^ ^ 



Let us now study the large-n behavior of the series ( [1.2|) . Using the expression with 
conjugate partitions (|2.8|) we find 



J2n 



if En 



'{n-\'^ + ti-l){n-X'^ + ti-2)...{n-\[ + ti-t + l)^ ^ 



n\) ^ J--^ \ {n + ti-l){n + ti-2)...in + ti-t+l] 



Since conjugate partitions are limited by the number n we can approximate the sum in 
( p.8|) for large n as an integi 
the first line Ai = p we get : 



for large n as an integral over the variables = Vi. Calling the number of boxes in 



2t-2 



The integrand is minimized by (1 — I'i)^*"^ and maximized by 1. In both cases one can 
symmetrize over the Vi to compute the integral explicitly and one finds : 

2 / r \2 



Therefore, the large-n behavior of the is fully controlled by the (cn/n!)^ prefactor and 



Z2n ~ exp 



2n ( ^ - 1 ) logn + 0(n) 



(2.11) 



An important conclusion is that for t > 2 the radius of convergence of (|2.1[) is infinite. 
Moreover approximating the sum over n by an integral we find 



Z exp constant x^-^ , (2-12) 



in agreement with (L6). This behavior is weU expected on physical grounds. Indeed the 
partition function reads also Z = exp — 5) where T = 1/ R is the temperature of the 
equivalent one-dimensional quantum system, and S and S are the impurity energy and 
entropy respectively. In the deep IR the impurity entropy converges to s ^ Inw and the 
energy, on dimensional grounds, scales as g'^ . The behavior (p. 121) is the analog of the 



"bulk term" in flows between bulk critical points. 

Although expressions (|2.7|) or (|2.8|) are in effect a solution of the problem, one can 
wonder about their practical use. Trying to evaluate the numerically, one finds that 
the series converges very slowly. For example, for t=3, evaluating a billion terms gives an 
accuracy of only about .1%. Fortunately, results are greatly improved by studying the free 
energy JF = —Tin 2, whose expansion we write 

oo 

^ = ^5^/2.0:2-, (2.13) 

n=0 

The /2n are of course given in terms of the Z2n- For example, /2 = —Z2, and = 
— Z4 + Z2/2. When evaluating the f2n numerically for n > 1, we find a much faster 
convergence. For t—3, we have 

/2 = -r(l/3)/r2(2/3) /s = 0.044223558 

/4 = 0.229454064 /lo = -0.022852208 (2.14) 

/e = -0.092261103 /12 = 0.012329254 

The Z2n can be extracted from these data, and are 

Z2 = r(l/3)/r2(2/3) Zs = 0.0618476490 

Z4 = 0.8378042270 Zio = 0.01021005440 (2.15) 

Za = 0.276783312 Z12 = 0.00131673987 

These coefficients are enough to get a good approximation of the properties all the way 
to the infrared (very large x) using Fade approximants. It is then preferable to consider 
the entropy S = which is bounded for x large. Keeping the coefficients through /12, 
one finds for instance that S{x = 0) — S{x = 00) ^ .57, in good agreement with the exact 
value Inv^^ .549306... 
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3. Other quantities of interest 

The previous calculation of the partition function is the simplest calculation which 
can be done using Jack symmetric functions. In this section we present several other 
calculations, and present a conjecture for the experimentally-measurable conductance. 

3.1. Twisted partition functions 

We have so far considered a periodic field (f) on the cylinder. We could also have 
winding modes such that 

(/)(a,r + i?) =(/>(a,r) + ^p, (3.1) 

where p is an integer. By splitting the field into classical and quantum parts we obtain an 
action similar to (|1.1|) but the interaction term is now 

dr cos B(h(n = D. -I- 2?: 

R J 



2g f dr cos \/3(j){a = 0, r) + 2n^T 
Jo R ■ 



Defining as before Z{p) = Z{g,p)/Z{g = 0,p) we find a perturbative expansion similar to 
( p. . 2|) with however each term in the sum multiplied by 

exp i2n^{Ti + . . . + Tn - t[ . . . - t'J . 

After change of variables we have the same expansion as but with 



7 . . 1 /ttM^^ dz[ \ iAiz)Aiz))HA{z')Aiz'))T f z^...Zr. Y 
^^-^P^- inlrfn[2^.z.2^nzJ U.^i^ ' - z'.z.p \z[...z'J ' ^'"'^ 

so in the Coulomb gas language there is now a magnetic charge located at the center of 
the circle. (We assume p is positive, otherwise just replace p by \p\.) We now use the fact 
that 

(zi . . . ZnYPx(z, a) = P\+p{z, a), (3.3) 
where X+p means the partition A where p columns of length n have been added. Therefore 

Z2n{p) = ^2 E ^>^^>^+P (3.4) 

!{A)<n 

where we use the relation Nx = A^a+p that follows immediately from the integral defining 
the norm N. The coefficients 6a read 



,,,A.-, + l + i(A^-z) 



For instance one has 



r(i + AOr(i + Ai+.) ^ sinf r(i-f) 

^^or2(i)r(i + Ai)r(i + Ai + ^) sin7r(i +^)r(i - i + ^)r(i - i - z)' 

(3.5) 

for X an arbitrary real number, and when a; = p is an integer, 

Up) = \i '\ r- (3.6) 

r(i - i +p)v{i - j-p) 

We can also compute the partition function when electrical neutrality is broken by 
some amount Q (assumed integer). Defining 

we have 

oo 
n=0 



where Zn^n+Q has formally the same expression as ( p.2| ) but there are n variables z and 



n-\- Q variables z' . By the same manipulations we obtain 



n!(n + Q)! ^ 

l{X)<n 



Observe that for n = we recover the well-known expression 

3.2. Correlation functions and the conductance 

The tool of Jack polynomials should allow the perturbative evaluation of correlation 
functions. However the calculation requires the knowledge of branching coefficients which 
are not known yet in general. We will illustrate this with an example. Consider the 
two-point function of the field (j) itself. By the same perturbative approach we find 

Wmr')), = {^(r)0(r')>„ + - ^ ' „V ' C„, (3.9) 

p=l ^ ' n=l 



where 



C = 

*-'np — 



fr (dz^dzTX (A(.)A(.))V^(A(.-)A(zO)V^ 
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The calculation can be easily done for the first term p = 1 using ( |3.3| ) . One finds following 
the same lines as in the previous paragraph 

Cm = 2{Z2{n) - Z2(0)) (3.10) 

For general p, we can still decompose the Coulomb-interaction term between positive and 
negative charges using Jack polynomials as in (|2.3| ). However, we cannot use the relation 
since we have the extra factors We can decompose [pT ,15| 



i "ixH. n(..j)[^i-^+(Ai-J + l)/a] 



[-1 

For our problem a = |^ = 1/^ and the only tableaux which contribute are the ones of the 
form Ai > A2, . . . > At, 1, . . . , 1. The Jack polynomials in (|3.11| ) multiply the ones in ( |2.3|) . 



We therefore end up with the problem of determining the coupling coefficients 

Unfortunately these coefficients are not yet known in general. 

The two-point function of the field ^ is especially useful in physical applications: in 
the 2d boundary problem the Kubo formula relates it to the conductance 0, while in 
dissipative quantum mechanics it is the mobility PJ^. The conductance at the Matsubara 
frequency = follows from 

= ^ /'''rfr'(</,(r)</>(r'))e^-(---'). (3.12) 



and the term is easily picked up in ( |3.10|) . The DC conductance then is obtained by 
analytically continuing ( p.l2|) to n = 0, leading to 



t t'^ n 



Using (|3.5| ) it is easy to perform the limit and one finds finally 



t r(i + i) 
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in agreement with integral done without the Jack functions 0. 

We cannot for the moment compute the two-point function to aU orders. However, 
we have the foUowing conjecture for G to all orders: 



^ 1 2a; ^2 , 



(3.14) 



p=0 



so the first few terms read then 
G 



\+'^ [Z'^x^ + 2(Z^ - Z2Z'^)x^ + 3(Z^ - Z'^Z^ - Za'i + Z'^Z^x'' + . . .) 



We define the continuation of Z2n{v) to real values of p by simple substitution in (|3.4|) , so 



2 = 



n! 



1^2 



A 

l{\)<n 



" 1 1 ' 

+ -(n - z +1)) - ^{K + 1 + -(n - i)) 



i=l 



(3.15) 



where ip{z) = r{z)' /r{z). We can investigate these numbers numerically. For t=3, ( |3.15| ) 
gives Z'^ = -2.64996, Z'^ = -2.351, and Z'^ = -0.964. The first agrees with ( |3l^ ), and the 
others are in very good agreement with the conductance calculated using the TBA in the 
next section. The conjecture ( p.l4|) appears to be a reasonable form of the Kubo formula 
0, but we have not succeeded in deriving it rigorously. 

To close this section we would like to remark that, although the theory of Jack sym- 
metric functions generally deals with rational values of t, we expect all formulas obtained 
above to hold for any t > 2 by naive substitution. 



4. Non-perturbative treatment 

The boundary problem (|1.3| ) is integrable |^^. One can thus find the exact 5" matrix 



for the quasiparticles of the problem [^^, and then use the thermodynamic Bethe ansatz 
2^,^ to compute the free energy and the conductance . In this section, we describe 



these results, and use them to derive a variety of functional equations. These functional 
relations give non-perturbative equations for the free energy, and allow one to derive simple 
but non-trivial relations among the coefficients Z2n- Another functional relation relates the 
conductance to the free energy, thus giving a new fluctuation-dissipation theorem for this 
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system. We also find expansion coefficients at t=3 numerically from the TBA as another 
check on our results. 

The starting point of the TBA is the quasiparticle description of a two-dimensional in- 
tegrable field theory. These quasiparticles scatter among themselves and off of the impurity 
with a known S matrix. At any value of t, the quasiparticle spectrum includes the soliton 
and antisoliton, which we label by + and — respectively. Moreover, at coupling t, there are 
t — 2 "breather" states in the spectrum. The energy and momentum of these left-moving 
massless particles are parametrized by rapidity variable 6, so E = —P = fij. g^p{—0), 
where fi+ = fx- = fx and fXj = 2(U sin(7r/2(t — 1)) for the breathers. We define the density 
of states Ur and the density of filled states pr for each quasiparticle species r. Periodic 
boundary conditions gives the functional of the pr- The free energy can be written 

in terms of these quantities; demanding it be at a minimum gives another set of relations 
which determine the densities. These relations are most conveniently written in terms of 
the functions er{0), which are defined by 

1 ^ Pr_ 

1 -|- e'^'- Ur 

Notice that if the particles are free, Pr/ur is the Fermi distribution function. However, for 
t 7^ 2, the particles are not free, and the are determined by the TBA equation 

e,{e) = ^e-' de'cprsiO - e') ln(l e-^«(^')), (4.1) 



s 

where the label s runs over breathers (1 . . .t — 2) and ±. The functions iprs for integer 
t are given in [^l| or [^. We will not need them, because here these equations can be 



written in a much simpler form 

{t-1 



OO 



where N^s is the incidence matrix of the following diagram 

12 s t-3 / 

O O 0---0 O t-2 



The dependences on the ratios //a//^ seem to have disappeared from (|^), but they appear 
as an asymptotic condition: the original equations (^T|) indicate that the solution must 
satisfy 



Cr —e-^ ase^ -OO. (4.3) 
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4.1. The partition function 

The impurity free energy is given in terms of e+: 



^TBA 



B 



T 



de 



t - 1 



ln(l + e^+(^)) 



2cos(7r/2(t- 1)) J 27rcosh[(t- 
where a = ln{T/TB) 0. The first piece is the non-analytic term (|1.4| ). Since the 



(4.4) 



same 



kernel appears in (f4.4| ) and (|4.2|) , and Z can be written in a simpler form for many of 
the t. Using the relation (|1.5| ) between JF and J^tba, we have for example 



= V>"i(«)/3 t = 3 

Z{a) = \{Y2{a)f'^ t = A 
Z{a) = ^Y^{a)/hY^{a) t = 5, 



(4.5) 



where we define Y^. = exp(e^). 

We derive simple functional relations for Z[a) by continuing it and the Yr{(y) into the 
complex ct-plane [^|. Using the simple identity 



lim 



A 



+ 



A 



cosh(A6' + i7r/2 — x) cosh(A6' — i7r/2 + x) 



Z(a + 7)2(« - 7) = ^(1 + > + («)), 



27r5(6'), 



we have 



where 7 = Z7r/2(t — 1). Similarly, the equations (|4.2|) yield 



(4.6) 



(4.7) 



Y+{e + 7)y+(^ - 7) = 1 + yt-2{e) 

Yt-2{9+-i)Yt-2{e--i) = {i+Yt-m){^+y+{d)f 

Ya{e + -i)Ya{e--f) = (i + y,+i(0))(i + y„_i(0)) 



(4.^ 



where a = 1 . . .t — 3, we define Yq = 0, and it follows from symmetry that Y^ = YL. 
These equations are applicable everywhere in the complex ^-plane, whereas the original 
TBA equations apply only in a strip \Im9\ < n/(t — l). The functional relations determine 
the functions Yr and Z uniquely once the asymptotic condition (|4.3|) is imposed. One can 



argue ||3^ that the functions Yr{a) (and Z{a)) have the periodicity Yr{a + t'j) = Yr{a), 



^ This can be derived from the kernels Ka of [p2| by using ( [4.2| ) along with the identity 
2cosh yKa = NabKb for a = 1 ... t — 2 and 2coshy(K-|-K_) = 2kt-2 + 1- 
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which imphes that they can be expanded in powers of (Tb/T)-2(*-i)/* for Tb/T smaU. 
This then gives the expansions (|2.1| ) and (|2.13D , because x oc (Tb/T)*^*"^)/*. 



Plugging the perturbative expansions into the functional relations (|4.7|) and (^3) give 
non-trivial relations among the coefficients determined in sect. 2 by the Jack-polynomial 
expansions. For t = 3, 4 these constraints can be written in a simple form by using ( [4.5|) . 
Doing a little algebra, we have 

3Z(a + iTi/2)Z{a - m/2)Z{a) = Z{a -f m/2) + Z{a) + Z{a - m/2) t = 3 

(4.9) 

4:Z{a + in/3)Z{a - i7r/3)Z{a) = Z{a + in/3) + 2Z^{a) + Z{a - z7r/3) t = 4. 

These relations have the nice feature that they have lost all trace of the quasiparticle index 
r. This is a strong hint that they can be derived directly, without having to do the full 
TBA analysis. One can also hope that there is a simple relation even for non-integer t 
(where the TBA analysis can get quite complicated); we make a conjecture in the next 



subsection. Even though (^4.9|) are stronger relations than (|4.7|) and (^]5) , this still probably 
isn't the end of the story: the conjecture for t = 4 amounts to 2Z{a)Z{a + 2in/3) = 
Z{a + iTv/3) + Z{a — iTv/3), which yields (|4.9|) but not the other way around. 
Plugging ( p.l|) into (|4.9|) , one finds 



exp -3 J2 /en^:'" j = ^6nX^" t = 3. 

\ n / n 

This means, for example, that for t=3, 2Zq = —6fe = 3Z2Z4 — [Z'lf' and Zyi = — 3/i2 + 
9(/6)^/2. Both agree with the Jack-polynomial expressions numerically evaluated in ( |2.14| ) 
and ( |2.15| ). In general, it means that that for t=3, the coefficients Z^n are given in 
terms of the lower coefficients. Similarly, for t=4 one finds that the coefficients Z^n are 
determined in terms of lower coefficients. For example, Z4 = (^2)^/3 = 7r/(3r^(3/4)) and 
9^8 = V^Z^Zq — Z|. It would certainly be interesting to have a direct proof (i.e. one 
depending only on the expressions (|2.7|) or ( |2.8| )) of these relations. They are certainly a 
hint of a much-deeper structure to the problem. 

Since we cannot determine all of the coefficients from the functional relations, a final 
check is to solve the TBA equations numerically and then fit the results to a power series. 
Doing a perturbative expansion of the non-perturbative solution, one obtains 

Ttba = tY,^iA-y) +-^- (4.10) 

n=0 ^ ^ 
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We can now match these results with those of the Jack-polynomial expansion, once com- 
parison of the first order has determined the ratio of the unknown constants k,' . We 
expect 

(4.11) 



/2n V /2 , 

We evaluate the full function /(T) to double- precision accuracy by solving the integral 
equations (|4.1|) , plugging this into the free energy (^^), and then fitting this to the series 
( [1.1 1|) at large T. For t=3, we find 



= -.4567084 /cg = 0.000422 
/C4 = 0.0224220 /cio = -0.00007 (4.12) 
/C6 = -0.002818 

To the accuracy of the TBA fit, we have excellent agreement. The scales x and Tg are 
therefore related for t = 3 by 

..Hm3H§(^V" (4.13) 

4-2. The conductance 



r(l/3) V T 



In this subsection we use the non-perturbative TBA to derive a remarkable fluctuation- 
dissipation relation of the conductance to the partition function. This allows us to obtain 
the value of infinitely-many coefficients in the perturbation expansion of G. It also allows 
us to conjecture a functional relation for Z for any rational value of t. 

The TBA gives the conductance as 

/■°° t-1 1 



oo 



2cosh2[(t-l)(6'-a)] l + l + (6')' 



where a = ln(T/TB) as before, and 1+ is given by the TBA equations (|4.2| ). By using the 
relation 



lim 



A2 



one finds 



cosh^(A6' -t- z7r/2 - x) cosh^(A6' - iTv/2 + x) , 

71 d 1 



G{a + 7) — G{a — 7) = — z 



t ~ldal + Y+{a) 
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where 7 = Z7r/(2(t — 1)). Using the relation ( [4.7|) , this gives 

G{a + -i) - G{a-^) = ~ij^^^^[Z-\a + ^)Z-\a~^)]. (4.15) 

This fluctuation-dissipation relation has lost all trace of the quasiparticles of the TBA: it 
is thus tempting to conjecture that it holds for all t, not just the integer values where the 
TBA analysis is valid. The perturbative expansion is Gpert = J2n d'^riX^"' ^ so we have for 
example 

27r TT 

92 = -^2^ cot y, 

in agreement with the perturbative calculations ( p.lO|) and (|3.13| ), which are valid for any 
t > 2. 

When t = p/q is rational, this gives many but not all of the g2n in terms of the Z2m 
(with m < n), because the terms on the left-hand side of ( [4.15| ) vanish when n is a multiple 
of p. For the physically-important value t=3, this gives all coefficients g6n+2 and g6n+4', 
for example 

2n 47r Stt 

92 = -^Z2 ^4 = ^^4 gs = ^{Zs-Z2Z,). 

These are in excellent agreement with a numerical calculation of the TBA conductance. 

Although the fact that some of the terms on the left-hand side of ( |4.15D vanish for 
rational t means we do not know how to relate these g2n to the Z2n, it does seem to imply 
a constraint on Z. We know that Z is an analytic function of x for t > 2, but we do not 
know that G is as well. (One way of checking this would be to check that the explicit 
perturbative expansions (|3.15| ) and ( p.7|) for G and Z obey the formula (|4.15| ).) If G is 



indeed analytic (so G=Gpert)^ then it requires that these terms on the right-hand-side also 
vanish, which means that Z should satisfy 

p 

J2 Z-\a -f 2j7)Z-i(a -f 2(j - 1)7) = p. (4.16) 
i=i 



For t=3 this is the relation already derived in (|4.9|) , but for t = 4 it is different. Putting 
it together with ([4.9| ) for t = 4, we find the simpler relation 2Z{a)Z[a + 2z7r/3) = Z{a + 
z7r/3) -|- Z[a — in /Z). This relation alone implies both ( [4.16| ) and ( [4.y| ) for t = 4. We can 



check the relationship ( [4.16| ) numerically. Plugging the expansion of Z into ( ^.16| ) gives 



the coefficients Z2p in terms of lower ones, which can be compared with the expression 
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(p77|). We have checked Zio for t = 5 and t = 5/2, and find that it is indeed satisfied. 
This leads us to conjecture that (|4.16 ) is true for aU t > 2 and rational. We also note 



that the relations (5]7) and ( f4.8| ) require that Z should obey an even more restrictive 
functional relationship. We have not succeeded in finding its general form, but one can 
always plug the perturbative expansions into (f4.7| ) and ( |4.8|) to derive more relations among 
the coefficients. 

To conclude this section, we recall first that the TBA analysis is usually made for t 
rational only. Moreover some of the results given above hold for t integer only. However 
some of the functional relations we have uncovered seem to make sense for any t. Observe 
also that the Jack expansion and TBA behave differently as t ^ 2. In the former case 
all integrals just blow up, while in the latter one gets finite results for the free energy, 
involving however logarithmic terms. Presumably, the TBA gives the regularized version 
of the Jack computations. 

5. Conclusion 

Using Jack polynomials and the thermodynamic Bethe ansatz, many properties of the 
Id log-sine gas can be computed exactly, some of which are of experimental significance. 
We hope that these methods can be used for other problems with potential applications, in 
particular for dissipative quantum mechanics. Moreover, multiple integrals similar to those 
we do using the Jack symmetric functions appear in many different kinds of computations, 
so we hope that these techniques are generally applicable. 

On the more formal side, it is exciting to have an example where two different areas of 
mathematical physics meet. By analogy, one might hope that these TBA techniques can be 
applied to other models, like the Calogero-Sutherland model, where Jack polynomials 
have been used recently. This overlap of techniques has led to intriguing relations between 
various quantities of Jack symmetric function theory. For example, when t=4 the series 
( p77| ) can be summed to give = 7r/(3r^(3/4)), but we have no direct proof, so we do not 
know if this is a fiuke or if a closed-form expression can be found for all t. The overlap has 
also led to several simple but powerful conjectures like ( |4.16| ) and ( |3.14| ). One can hope 



that this is evidence of a more complete mathematical structure behind the scene. 
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